Impacts of preference and geography on epidemic spreading 
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We investigate the standard susceptible-infected-susceptible model on a random network to study the effects 
of preference and geography on diseases spreading. The network grows by introducing one random node with 
m links on a Euclidean space at unit time. The probability of a new node i linking to a node j with degree kj at 
distance dij from node i is proportional to /df^, where A and B are positive constants governing preferential 
attachment and the cost of the node-node distance. In the case of yl = 0, we recover the usual epidemic behavior 
with a critical threshold below which diseases eventually die out. Whereas for _B = 0, the critical behavior is 
absent only in the condition A — 1. While both ingredients are proposed simultaneously, the network becomes 
robust to infection for larger A and smaller B. 

PACS numbers: 89.75.Hc, 87.23.Ge, 05.70.Ln, 87.19.Xx 



The classical mathematical approach to diseases spreading 
either ignores the population structure or treats populations as 
distributed in a regular medium iH, 0] . However, it has been 
suggested recently that many social, biological, and commu- 
nication systems possess two universal characters, the small- 
world effect 1 3] and the scale-free property |4], which can be 
described by complex networks whose nodes represent indi- 
viduals and links represent interactions among them |5, 6]. In 
view of the wide occurrence of complex networks in nature, 
it is important to study the effects of topological structures 
on the dynamics of epidemic spreading. Pioneering works 
QBHIlSIllll have given some valuable insights: for homo- 
geneous networks (e.g., exponential networks), there are crit- 
ical thresholds of the spreading rate below which infectious 
diseases will eventually die out; on the contrary, even infec- 
tions with low spreading rates will prevail over the entire pop- 
ulation in heterogeneous networks (e.g., scale-free networks). 
This radically changes many conclusions drawn from clas- 
sic epidemic modelling. Furthermore, it has been observed 
that the heterogeneity of a population network in which the 
disease spreads may have noticeable effects on the evolution 
of the epidemic as well as the coiTesponding immunization 
strategies lIllQlHEIll. 

For many real networks, however, individuals are embed- 
ded in a Euclidean space and the interactions among them 
usually depend on their sgatial distances and take place among 
their nearest neighbors [ O, um flStl . For instance, the number 
of long-ranged links and the number of edges connected to a 
single node are limited by the spatial embedding, particularly 
in planar networks. Preferential attachment is weakened by 
geographical embedding |18]. Also, people have proved that 
the characteristic distance plays a crucial role in the dynamics 
taking place on those networks ll20l l2ll l22i l23ll . Thus, it is 
natural to study associated influences of preference and geog- 
raphy on epidemic spreading. But up to now only a few of 
works address this problem, e.g., modeling transmission as a 
function of geographical distance l24il25ll availably capturing 
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FIG. 1: Degree distributions of the generated networks with mo = 
m = 3 in two special cases, ^ = (a) or B = (b). The size of the 
network is Ai' = 10^. 



the dynamics of diseases in wild and domesticated animals 
MM- 

In this paper, we study the standard susceptible-infected- 
susceptible (SIS) model on a growing network in Euclidean 
space. On a vertical plane, the growth of the network depends 
jointly on two mechanisms, preference and geography. The 
placement of links is driven by competition between preferen- 
tial attachment and distance dependence. In the case that the 
network grows with geographical constraint, we recover the 
usual epidemic behavior with a critical threshold below which 
diseases will eventually die out. While the network is totally 
governed by preferential attachment, the epidemic behavior 
depends on the preferential exponent. When both factors are 
considered simultaneously, it becomes difficult for epidemic 
spreads as the preference has an overwhelming majority than 
the geography. 

Specifically in two-dimensional x — y plane, we consider 
a square of unit size and with periodic boundary conditions. 
To construct a network of N nodes, let {xi,X2, ■ ■ ■ ,xn) 
and {yi,y2, ■ ■ ■ jVn) be the 27V independent random vari- 
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TABLE I: The average path length L and clustering coefficient C of 
the generated networks with size A'^ = 10000 in the case of A = 
for different values of B. 
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L 


C 


0.5 


5.13(1) 


0.00(1) 


1 


5.13(5) 


0.00(2) 


2 


5.40(1) 


0.03(9) 


3 


6.41(7) 


0.19(2) 


4 


7.11(2) 


0.31(0) 



ables identically and uniformly distributed within the inter- 
val {0, 1}. A specific set of values of the random variables 
{{xi,yi), {x2,y2), {xN,yN)} is chosen to represent the 
coordinates of the N randomly distributed nodes Q. The 
network starts with mo nodes and then the other nodes with 
TO links are added one by one at each time step according to 
their serial numbers i = toq + 1 to A^. Following ideas pro- 
posed by Yook et al. flS*], the probability that a new node i 
links to a old node j with kj links at distance dij from node i 
is 

l[^h,d,j)-%, (1) 

"'ij 

where A and B are positive constants, governing preferen- 
tial attachment and the cost of the node-node distance. We 
note following interesting features, (i) In the case of A — 0, 
the network is geographically grown with an exponential dis- 
tribution of nodes' degree (see Fig. In the limit of 
B ^ (X, only the smallest value of d coii'esponding to the 
nearest node will contribute with probability 1 llisll . (ii) In the 
case of i? = 0, the network reduces to the Barabasi-Albert 
(BA) graph only for A — 1. In the region < A < 1, 
the nodes' degree distribution is stretched exponential. For 
A > 1, a finite number of nodes connect to nearly all other 
nodes Q. That is illustrated by Fig. [Itb). 

To estimate the effect of the network's topology on epi- 
demic dynamics, we will investigate the standard SIS model 
This model relies on a coarse-grained description of indi- 
viduals in the population. Namely, each node of the network 
represents an individual and each link is a connection along 
that the infection can spread to other individuals. The individ- 
uals can only exist in two states, susceptible and infected. At 
each time step, each susceptible node is infected with proba- 
bility u if it is connected to one or more infected nodes. At the 
same time, the infected nodes become susceptible again with 
probability S, defining an effective spreading rate A = ly/S. 
We can set 5 = 1 without lack of generality, since it only af- 
fects the definition of the time scale of the virus propagation. 
Individuals run stochastically through the cycle susceptible 
infected susceptible. 

Let us focus on the case A = Q first, i.e., the network 
grows with geographical constraint. According to Eq. ([T]|, 
the preferential attachment is excluded. Since all nodes are 
uniformly distributed in the square, the only effect of the fac- 
tor B is determining the average path length of the network 
while the degree distribution of nodes has the same behav- 




FIG. 2: Density of infected nodes p as a function of time t in the 
network with A = and B — 1, 2, 3, respectively. The spreading 
rate is A = 0.15. Simulations were performed on networks with 
N = 10000 and m = mo = 3. 



0.4 






1 ' 1 ' 
«'» 

a' 










0.3 






0.2 






/ : 


0.1 




1/ 


-□-B = 0.5 - 
-o-B=1 


0.0 






-o-B = 2 
-H-B=4 
I.I. 




0.0 


0.2 


0.4 0.6 



X 



FIG. 3: Density of infected nodes p as a function of A in the network 
with A = Q and B = 0.5, 1, 2, 4, respectively. 



ior. For small B, the role of node-node distance is weak and 
old nodes are linked with approximate randomness. When B 
becomes large, the geographical influence is strong and only 
nodes around the new one will be connected with large pos- 
sibility, hence the local clustering. This feature is reflected in 
Tab. U that is, the average path length L and the clustering 
coefficient C get larger with the increase of B. 

Figure |2] shows the evolution of the infected nodes density 
as a function of time for epidemics with A = 0.15. We start 
from a single infected node of the network, and iterate the 
rules of the SIS model with parallel updating. Each curve 
represents the average over 10 different starting configura- 
tions, performed on 10 different realizations of the random 
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networks. We clearly notice a great influence of the geogra- 
phy on the spreading velocity of diseases, namely, the smaller 
the parameter B is, the more fast the infection propagates. For 
finite B, since the geography does not change the network's 
connectivity distribution (see Fig. [11^))^ the evolution of node 
i's degree can be written as 



dki(t) m 
~~dt T' 



(2) 



with the initial condition ki{i) = m. One can easily write the 
solution 



hit) = m(ln- + 1), 
and accordingly obtain the degree distribution 



P{k) = \j25{Ht)-k)^^ 



(3) 



(4) 



i=l 



In complex networks, the basic reproductive number takes the 
form, i?o ~ [13] ■ Different from the classical result, 

it defines an epidemic threshold Ac — (k) / (fc^). Combining 
Eq. (|4|l, we have 



J m 

J m 



2 

5m 



(5) 



In Fig. [3] we plot the steady density of infected nodes p as a 
function of the spreading rate A for the case of A = Q, which is 
the time average of the fraction of infected individuals reached 
after an initial transient regime. Simulations were computed 
over 50 different starting configurations, performed on 50 dif- 
ferent realizations of the networks. The size of networks is 
N — 10^. As shown in Fig. [3] all curves display the same 
behavior and the SIS model exhibits an epidemic threshold, 
Ac = 0.12(2), despite the variety of B, which is in agreement 
with the analytical prediction, Ac = 2/(5 x 3) = 0.133 (Eq. 
©) H^. That is different from the results gained by Santos et 
al., who also studied the SIS model on a homogeneous small- 
world network f3^, the critical value of A changes smoothly 
as one varies rewiring probability without changing the degree 
distribution. In our model, the increase of B gives rise to the 
average path length and the clustering coefficient. Whereas in 
isoll . those network features reduce with the increase of the 
rewiring probability. 

Next, we will investigate the dynamics in the case of _B = 0, 
i.e., the network grows following preferential attachment. For 
the linear preference, A — 1, the generated network reduces 
to the exact BA graph. As shown in Fig. 21 the behavior of the 
density of infected nodes follows the property, p ^ e^^^™^ 
inil . which implies the surprising absence of any epidemic 
threshold in the model, i.e., Ac = 0. As to the nonlinear case, 
we first consider the region A < 1 and go down a little from 
A = 1. In this event, highly connected nodes become less 
attractive for attachment compared with the linear preference. 
The resulting degree distribution is of the form |28], P{k) ^ 
exp{— iik^^^ / {1 — A)), where /i is a positive constant 
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FIG. 4: Density of infected nodes p vs. A (a) and 1/A (b) in the net- 
work with B = 0. The preference exponents are A = 0.5 (closed 
squares), 1 (open circles), 1.5 (closed diamonds), and 2 (open trian- 
gles), respectively. 



0.006875 
0.01375 
0.02063 
0.02750 
0.03438 
0.04125 
0.04813 
O.OSOOO 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

A 



FIG. 5: Density of infected nodes pin A~ B plane for the spreading 
rate A = 0.15 (a) and 0.5 (b), respectively. 
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depends on A, fi = Y^'^=2 11^=2(1 + m/*"^)"^- I" the region 
A > 1, the attractiveness of the old highly connected nodes 
increases which results in a small number of nodes that get 
all connections in the network. As shown in Fig. H] there 
are epidemic thresholds of the SIS model for A = 0.5, 1.5, 
and 2, respectively. Furthermore, the threshold increases as A 
becomes larger 

Finally, we plot the prevalence p 'm A ~ B plane in Fig. |5] 
to shown the influence of the competition of two ingredients 
on infections. For A = 0.15, there exists a set of peaks in 
Fig. |5ja). Namely, given the value of B, as A is increased 
initially, the density of infected nodes p increases gradually 
and reaches a maximum for some value of A, and then de- 
creases rapidly to as A is increased further. Furthermore, 
the contour planes of p takes an excursion to right with the 
increase of B. We argue that the following factor should be 
taken into account to understand this performance. According 
to Eq. ([Til, as B increases, the new nodes are preferential to 
connect their nearest neighbors, which results in the decrease 
of the number of long-ranged links. The network becomes 
more local clustering and robust to epidemic spreading. To 
keep the same prevalence, the effect of node's degree should 
be strengthened, i.e., increases A. Thus the contour planes 
lean to right as B becomes larger For A — 0.5, the prevalence 



p displays a different behavior As shown in Fig. |3b), p de- 
creases monotonically as A gets larger, and if B decreases at 
the same time, the network becomes robust to diseases. 

To summarize, we have studied the SIS model on a random 
network. On an Euclidean plane, the network grows depend- 
ing jointly on the preference and the geography. The former 
indicates the attractiveness of highly connected nodes and the 
later denotes the geographical constraint. It is found that both 
factors have great influences on the infection. For the net- 
work growing with the geographically constraint {A = 0), we 
recover the usual epidemic behavior with a critical threshold 
below which diseases will eventually die out. While the net- 
work is purely governed by preferential attachment {B — 0), 
the epidemic behavior depends on the preferential exponent 
and the critical phenomenon is absent only in the condition 
A — 1. When both factors are present simultaneously, the 
network becomes robust to diseases as the preference has an 
overwhelming majority than the geography. In real world, 
agents located on different positions according to the compe- 
tition between the preference and the geography. The above 
description of the spreading dynamics might contribute to un- 
derstanding realistic epidemics. 
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